Identification of the novel null allele HLA-DPB1*1449:01N by next-generation sequencing.

The novel HLA-DPB1*1449:01N allele differs from HLA-DPB1*16:01:01:01 by a nucleotide at codon 92 in exon 2.

The complex HLA-E-nonapeptide in Behçet disease.

Introduction: The knowledge of the aetiology of Behçet disease (BD), an immune-mediated vasculitis, is limited. HLA-B, mainly HLA-B51, and HLA-A molecules are associated with disease, but the ultimate cause of this association remains obscure. There is evidence that NK cells participate in the etiopathology of BD. NK cells have activator and inhibitor surface receptors, like the KIR and the NKG2 families. Classical HLA-class I molecules (A, B and C) are keys in the activity control of the NK because they are KIR ligands. Most NKG2 receptors bind HLA-E, which presents only nonapeptides derived from the signal peptide of other class-I molecules. Objective: This study investigates the contribution of the pair HLA-E and ligand, nonapeptide derived from the 3-11 sequence of the signal peptides of class I classical molecules, to the susceptibility to BD. Methods: We analyzed the frequency of the HLA-derivated nonapeptide forms in 466 BD patients and 444 controls and an HLA-E functional dimorphism in a subgroup of patients and controls. Results: In B51 negative patients, the frequency of VMAPRTLLL was lower (70.4% versus 80.0% in controls; P=0.006, Pc=0.04, OR=0.60, 95%CI 0.41-0.86), and the frequency of VMAPRTLVL was higher (81.6% versus 71.4% in controls; P=0.004, Pc=0.03, OR=1.78, 95%CI 1.20-2.63). In homozygosity, VMAPRTLLL is protective, and VMAPRTLVL confers risk. The heterozygous condition is neutral. There were no significant differences in the distribution of the HLA-E dimorphism. Discussion: Our results explain the association of BD with diverse HLA-A molecules, reinforce the hypothesis of the involvement of the NK cells in the disease and do not suggest a significant contribution of the HLA-E polymorphism to disease susceptibility.

Monitoring of kinetics and exhaustion markers of circulating CAR-T cells as early predictive factors in patients with B-cell malignancies.

Purpose: CAR-T cell therapy has proven to be a disruptive treatment in the hematology field, however, less than 50% of patients maintain long-term response and early predictors of outcome are still inconsistently defined. Here, we aimed to optimize the detection of CD19 CAR-T cells in blood and to identify phenotypic features as early biomarkers associated with toxicity and outcomes. Experimental design: In this study, monitoring by flow cytometry and digital PCR (dPCR), and immunophenotypic characterization of circulating CAR-T cells from 48 patients treated with Tisa-cel or Axi-cel was performed. Results: Validation of the flow cytometry reagent for the detection of CAR-T cells in blood revealed CD19 protein conjugated with streptavidin as the optimal detection method. Kinetics of CAR-T cell expansion in blood confirmed median day of peak expansion at seven days post-infusion by both flow cytometry and digital PCR. Circulating CAR-T cells showed an activated, proliferative, and exhausted phenotype at the time of peak expansion. Patients with increased expansion showed more severe CRS and ICANs. Immunophenotypic characterization of CAR-T cells at the peak expansion identified the increased expression of co-inhibitory molecules PD1 and LAG3 and reduced levels of the cytotoxicity marker CD107a as predictors of a better long-term disease control. Conclusions: These data show the importance of CAR-T cells in vivo monitoring and identify the expression of PD1LAG3 and CD107a as early biomarkers of long-term disease control after CAR-T cell therapy.

Factors associated with the humoral response after three doses of COVID-19 vaccination in kidney transplant recipients.

Introduction: Kidney transplant recipients showed a weak humoral response to the mRNA COVID-19 vaccine despite receiving three cumulative doses of the vaccine. New approaches are still needed to raise protective immunity conferred by the vaccine administration within this group of high-risk patients. Methods: To analyze the humoral response and identify any predictive factors within these patients, we designed a prospective monocentric longitudinal study of Kidney transplant recipients (KTR) who received three doses of mRNA-1273 COVID-19 vaccine. Specific antibody levels were measured by chemiluminescence. Parameters related to clinical status such as kidney function, immunosuppressive therapy, inflammatory status and thymic function were analyzed as potential predictors of the humoral response. Results: Seventy-four KTR and sixteen healthy controls were included. One month after the administration of the third dose of the COVID-19 vaccine, 64.8% of KTR showed a positive humoral response. As predictive factors of seroconversion and specific antibody titer, we found that immunosuppressive therapy, worse kidney function, higher inflammatory status and age were related to a lower response in KTR while immune cell counts, thymosin-a1 plasma concentration and thymic output were related to a higher humoral response. Furthermore, baseline thymosin-a1 concentration was independently associated with the seroconversion after three vaccine doses. Discussion: In addition to the immunosuppression therapy, condition of kidney function and age before vaccination, specific immune factors could also be relevant in light of optimization of the COVID-19 vaccination protocol in KTR. Therefore, thymosin-a1, an immunomodulatory hormone, deserves further research as a potential adjuvant for the next vaccine boosters.

Higher plasma levels of thymosin-α1 are associated with a lower waning of humoral response after COVID-19 vaccination: an eight months follow-up study in a nursing home.

BACKGROUND: Older people achieve lower levels of antibody titers than younger populations after Covid-19 vaccination and show a marked waning humoral immunity over time, likely due to the senescence of the immune system. Nevertheless, age-related predictive factors of the waning humoral immune response to the vaccine have been scarcely explored. In a cohort of residents and healthcare workers from a nursing home that had received two doses of the BNT162b2 vaccine, we measured specific anti-S antibodies one (T1), four (T4), and eight (T8) months after receiving the second dose. Thymic-related functional markers, including thymic output, relative telomere length, and plasma thymosin-α1 levels, as well as immune cellular subsets, and biochemical and inflammatory biomarkers, were determined at T1, and tested for their associations with the magnitude of the vaccine response (T1) and the durability of such response both, at the short- (T1-T4) and the long-term (T1-T8). We aimed to identify age-related factors potentially associated with the magnitude and persistence of specific anti-S immunoglobulin G (IgG)-antibodies after COVID-19 vaccination in older people. RESULTS: Participants (100% men, n = 98), were subdivided into three groups: young (< 50 years-old), middle-age (50-65 years-old), and older (≥65 years-old). Older participants achieved lower antibody titers at T1 and experienced higher decreases in both the short- and long-term. In the entire cohort, while the magnitude of the initial response was mainly associated with the levels of homocysteine [ß (95% CI); - 0.155 (- 0.241 to - 0.068); p = 0.001], the durability of such response at both, the short-term and the long-term were predicted by the levels of thymosin-α1 [- 0.168 (- 0.305 to - 0.031); p = 0.017, and - 0.123 (- 0.212 to - 0.034); p = 0.008, respectively]. CONCLUSIONS: Higher plasma levels of thymosin-α1 were associated with a lower waning of anti-S IgG antibodies along the time. Our results suggest that plasma levels of thymosin-α1 could be used as a biomarker for predicting the durability of the responses after COVID-19 vaccination, possibly allowing to personalize the administration of vaccine boosters.

Sex-specific analysis in Behçet's disease reveals higher genetic risk in male patients.

OBJECTIVES: Behçet's disease tends to be more severe in men than women. This study was undertaken to investigate sex-specific genetic effects in Behçet's disease. METHODS: A total of 1762 male and 1216 female patients with Behçet's disease from six diverse populations were studied, with the majority of patients of Turkish origin. Genotyping was performed using an Infinium ImmunoArray-24 BeadChip, or extracted from available genotyping data. Following imputation and extensive quality control measures, genome-wide association analysis was performed comparing male to female patients in the Turkish cohort, followed by a meta-analysis of significant results in all six populations. In addition, a weighted genetic risk score for Behçet's disease was calculated and compared between male and female patients. RESULTS: Genetic association analysis comparing male to female patients with Behçet's disease from Turkey revealed an association with male sex in HLA-B/MICA within the HLA region with a GWAS level of significance (rs2848712, OR = 1.46, P = 1.22 × 10-8). Meta-analysis of the effect in rs2848712 across six populations confirmed these results. Genetic risk score for Behçet's disease was significantly higher in male compared to female patients from Turkey. Higher genetic risk for Behçet's disease was observed in male patients in HLA-B/MICA (rs116799036, OR = 1.45, P = 1.95 × 10-8), HLA-C (rs12525170, OR = 1.46, P = 5.66 × 10-7), and KLRC4 (rs2617170, OR = 1.20, P = 0.019). In contrast, IFNGR1 (rs4896243, OR = 0.86, P = 0.011) was shown to confer higher genetic risk in female patients. CONCLUSIONS: Male patients with Behçet's disease are characterized by higher genetic risk compared to female patients. This genetic difference, primarily derived from our Turkish cohort, is largely explained by risk within the HLA region. These data suggest that genetic factors might contribute to differences in disease presentation between men and women with Behçet's disease.

Protein and functional isoform levels and genetic variants of the BAFF and APRIL pathway components in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is the prototype of an autoimmune disease. Belimumab, a monoclonal antibody targets BAFF, is the only biologic approved for SLE and active lupus nephritis. BAFF is a cytokine with a key-regulatory role in the B cell homeostasis, which acts by binding to three receptors: BAFF-R, TACI and BCMA. TACI and BCMA also bind APRIL. Many studies reported elevated soluble BAFF and APRIL levels in the sera of SLE patients, but other questions about the role of this system in the disease remain open. The study aimed to investigate the utility of the cytokine levels in serum and urine as biomarkers, the role of non-functional isoforms, and the association of gene variants with the disease. This case-control study includes a cohort (women, 18-60 years old) of 100 patients (48% with nephritis) and 100 healthy controls. We used ELISA assays to measure the cytokine concentrations in serum (sBAFF and sAPRIL) and urine (uBAFF and uAPRIL); TaqMan Gene Expression Assays to quantify the relative mRNA expression of ΔBAFF, ßAPRIL, and εAPRIL, and next-generation sequencing to genotype the cytokine (TNFSF13 and TNFSF13B) and receptor (TNFRSF13B, TNFRSF17 and TNFRSF13C) genes. The statistical tests used were: Kruskal-Wallis (qualitative variables), the Spearman Rho coefficient (correlations), the Chi-square and SKAT (association of common and rare genetic variants, respectively). As expected, sBAFF and sAPRIL levels were higher in patients than in controls (p ≤ 0.001) but found differences between patient subgroups. sBAFF and sAPRIL significantly correlated only in patients with nephritis (rs = 0.67, p ≤ 0.001) and ßAPRIL levels were lower in patients with nephritis (p = 0.04), and ΔBAFF levels were lower in patients with dsDNA antibodies (p = 0.04). Rare variants of TNFSF13 and TNFRSF13B and TNFSF13 p.Gly67Arg and TNFRSF13B p.Val220Ala were associated with SLE. Our study supports differences among SLE patient subgroups with diverse clinical features in the BAFF/APRIL pathway. In addition, it suggests the involvement of genetic variants in the susceptibility to the disease.

HLA-E restricted cytomegalovirus UL40 peptide polymorphism may represent a risk factor following congenital infection.

BACKGROUND: Congenital cytomegalovirus immunopathogenesis is largely unknown and multifactorial due to the complex interactions between viral, maternal, placental, and child factors. Polymorphisms in the HLA-E binding UL4015-23 peptide mimics HLA-E complexed peptides from certain HLA-A, -B, -C and -G alleles, which regulate the cellular immune response driven by natural killer-cells (NK) and CD8 + T cells. The aim of this study was to compare UL4015-23 peptides distribution in congenital CMV and the counterpart HLA Class I peptides in a healthy cohort to investigate risk factors and markers for cCMV disease. In this 10-year retrospective study, the UL40 gene was directly sequenced from 242 clinical samples from 199 cases of congenital CMV (166 children and 33 pregnant or breast feeding women). Distribution of HLA-E binding UL4015-23 peptides was analyzed and compared to those of HLA Class I observed in a cohort of 444 healthy individuals. RESULTS: Nineteen different HLA-E binding UL4015-23 peptides were found. Three of them (VMAPRTLIL, VMAPRTLLL, VMAPRTLVL) were found in 88.3% of UL40 and 100% of HLA Class I of healthy individuals. In contrast, 15 of them (10.7%) were not found in HLA Class I. The VMAPRTLFL peptide was found in 1% of UL40 and all HLA-G alleles. Significant differences in peptide (VMAPRTLIL, VMAPRTLLL, VMAPRTLVL, other UL4015-23 peptides, other HLA Class I peptides) distribution between UL4015-23 from congenital CMV and HLA-A, -B, -C and -G from healthy individuals were found. CONCLUSIONS: Our findings suggest that a mismatch between UL4015-23 peptides and HLA Class I peptides between children and mothers might play a role in congenital CMV disease, and it may account for differences in outcome, morbidity and sequelae.

Variants in autoinflammatory diseases-related genes in a family with cold-induced autoinflammatory syndrome.

OBJECTIVES: To characterise what immunogenetic alterations are present in a Spanish family having several members with a familial cold-induced autoinflammatory syndrome (FCAS), a kind of autoinflammatory disease (AID). METHODS: We present the case of two sisters (cases 1 and 2) with a similar clinical picture since their childhood. The symptoms start after exposure to cold and consist of recurrent fever, papules or urticaria, and oedema in hands and fingers. The mother had similar symptomatology as her daughters, which remitted after her first pregnancy, whereas the father is healthy. Patients and their parents were genotyped in a panel of 14 candidate genes using Next-Generation Sequencing (NGS). Real-time PCR was used to quantify IL1ß mRNA levels from LPS-stimulated monocytes. ELISA was used to measure the IL1ß and IL18 concentrations in supernatants and sCD25 levels in sera. IL1ß, IL4, IL6, IL8, IL10, IL17A, IL18 and TNF-α serum levels were assessed using xMAP® Technology. RESULTS: All the genetic variants found in this family are benign with two exceptions: NLRC4 p.Leu339Pro (present in both cases and their mother) and PSTPIP1 p.Gln219His (present in Case 1 and her father). The monocytes stimulated of the individuals with the NLRC4 variant produce higher levels of IL1ß (protein and mRNA). Levels of TNF-α, IL4, and IL6 were higher in Case 1 than in the age-matched controls. CONCLUSIONS: The familial segregation and the clinical picture compatible with FCAS suggest that NLRC4 p.Leu339Pro causes the AIDs syndrome diagnosed in several family members.

Genetic Association of a Gain-of-Function IFNGR1 Polymorphism and the Intergenic Region LNCAROD/DKK1 With Behçet's Disease.

OBJECTIVE: Behçet's disease is a complex systemic inflammatory vasculitis of incompletely understood etiology. This study was undertaken to investigate genetic associations with Behçet's disease in a diverse multiethnic population. METHODS: A total of 9,444 patients and controls from 7 different populations were included in this study. Genotyping was performed using an Infinium ImmunoArray-24 v.1.0 or v.2.0 BeadChip. Analysis of expression data from stimulated monocytes, and epigenetic and chromatin interaction analyses were performed. RESULTS: We identified 2 novel genetic susceptibility loci for Behçet's disease, including a risk locus in IFNGR1 (rs4896243) (odds ratio [OR] 1.25; P = 2.42 × 10-9 ) and within the intergenic region LNCAROD/DKK1 (rs1660760) (OR 0.78; P = 2.75 × 10-8 ). The risk variants in IFNGR1 significantly increased IFNGR1 messenger RNA expression in lipopolysaccharide-stimulated monocytes. In addition, our results replicated the association (P < 5 × 10-8 ) of 6 previously identified susceptibility loci in Behçet's disease: IL10, IL23R, IL12A-AS1, CCR3, ADO, and LACC1, reinforcing the notion that these loci are strong genetic factors in Behçet's disease shared across ancestries. We also identified >30 genetic susceptibility loci with a suggestive level of association (P < 5 × 10-5 ), which will require replication. Finally, functional annotation of genetic susceptibility loci in Behçet's disease revealed their possible regulatory roles and suggested potential causal genes and molecular mechanisms that could be further investigated. CONCLUSION: We performed the largest genetic association study in Behçet's disease to date. Our findings reveal novel putative functional variants associated with the disease and replicate and extend the genetic associations in other loci across multiple ancestries.

Brief Report: Toll-like Receptor 9-1635A/G Polymorphism Is Associated With HIV-1 Rebound After Four Weeks of Interruption of Antiretroviral Therapy.

OBJECTIVES: This study aims to analyze the association of the presence of common polymorphisms [single nucleotide polymorphisms (SNPs)] on Toll-like receptors (TLRs), such as TLR9-1635A/G, TLR2-1892A/C, TLR2-2258G/A, TLR4-899A/G, and TLR4-1196C/T, with the viral rebound after stopping antiretroviral treatment (ART). CCR5-Δ32 deletion and HLA-A/HLA-B alleles were also analyzed. DESIGN: Interruption of ART may be required to investigate the outcome of strategies aimed to achieve drug-free HIV remission or cure. However, interruption of ART is currently not indicated. This was a retrospective longitudinal study that included 57 long-term suppressed HIV-1-infected individuals. METHODS: TLR SNPs were detected by real-time polymerase chain reaction (PCR). CCR5-Δ32 was analyzed by conventional PCR and HLA-A and HLA-B alleles by PCR-SSOP Luminex. RESULTS: HIV-1 RNA rebound at week 4 after treatment interruption positively correlated with pre-ART HIV-1 load (P = 0.025). The TLR9-1635AA genotype was independently associated with a higher HIV-1 rebound compared with those with AG + GG genotype (multivariate stepwise regression analysis, P = 0.004). Women had lower HIV-1 RNA load both at rebound and during the 72 weeks of follow-up compared with men (P < 0.05 at all time-points), whereas CD4 nadir and CD4 count set-point were similar according to sex. The pre-ART viral load was independently associated with the viral set-point (P = 0.001), whereas the presence of the HLA-A01 allele (P = 0.027) and the CD4 nadir (P = 0.001) were associated with the CD4 count set-point. CONCLUSIONS: The association of the TLR9-1635AA genotype with a higher HIV-1 rebound suggests that this SNP may affect the results from strategies requiring interruption of ART aimed to cure HIV-1 infection.

Lower frequency of anti-citrullinated protein antibodies among early arthritis patients with high body mass index.

OBJECTIVES: To investigate the role of body mass index (BMI) in the phenotypic and genotypic characteristics of early arthritis patients. METHODS: We analysed the clinical and laboratory parameters from the baseline visit of patients (670 patients [78.51% women]) included in the PEARL study. The WHO definition for low weight, normal weight, overweight and obesity (BMI <18.5, 18.5-25, 25-30 or ≥30 kg/m2, respectively) was applied. Anticitrullinated protein antibodies (ACPA) were studied by ELISA and HLA-DRB1* were genotyped by sequence speci c oligonucleotide probes. The relationship between BMI classification and other variables was analysed using Kruskall-Wallis, Anova and Chi-Square tests. Then multivariate logistic regression was performed to establish the role of BMI in ACPA positivity and ordered logistic regression to establish its relationship with ACPA level. RESULTS: Among the patients studied, 255 (38.06%) were considered overweight and 136 (20.3%) obese. High BMI patients had significantly more pain perception and disability than normal weight patients, whereas no clear differences in disease activity were observed between high BMI and normal weight patients. ACPA positivity was significantly less frequent in overweight and obese patients compared to normal BMI patients. This information was confirmed by adjusting for smoking habit and the presence of shared epitope. CONCLUSIONS: Our data support the theory that high BMI patients suffer more frequently from ACPA-negative RA. Nevertheless, although no disease activity differences were observed, these patients showed higher pain and disability scores since the beginning of disease.

The role of a functional variant of TYK2 in vasculitides and infections.

OBJECTIVES: The TYK2 gene encodes a tyrosin kinase which is involved in multiple immune functions. A functional variant of this gene has been identified to play a protective role in multiple autoimmune diseases. The goal of this study was to evaluate the involvement of this variant of TYK2 in vasculitides [giant cell arteritis (GCA), ANCA-associated vasculitis (AAV) and IgA vasculitis (IgAV)] and viral infections [hepatitis C virus (HCV) and human immunodeficiency virus type I (HIV-1)]. METHODS: The study sample was composed of 13,745 European individuals. The genotyping was performed by Immunochip and TaqMan 5' allele discrimination assays and the allele frequencies were compared using PLINK. RESULTS: Although the results obtained did not reach the genome-wide level of significance, p-values at nominal significance were observed, suggesting that the TYK2 variant provides protection against two vasculitides: GCA (p=5.94E-3; OR (95%CI) = 0.56 (0.37-0.85) and AAV (p=6.79E-3; OR (95%CI) = 0.65 (0.47-0.89). However, this variant was not found to be associated with IgAV. No evidence was gained that the TYK2 variant confers susceptibility to HCV and HIV-1 infection. CONCLUSIONS: This is the first study to propose the association between the TYK2 and both GCA and AAV. Our findings also suggest that TYK2 does not play a relevant role in IgAV or in susceptibility to HCV and HVI-1.

Characterization of a novel HLA-DRB1*03 allele, DRB1*03:171.

DRB1*03:171 was described in a Caucasian Spanish family.

Association of Functional Polymorphisms of KIR3DL1/DS1 With Behçet's Disease.

Behçet's disease (BD) is an immune-mediated vasculitis related to imbalances between the innate and adaptive immune response. Infectious agents or environmental factors may trigger the disease in genetically predisposed individuals. HLA-B51 is the genetic factor stronger associated with the disease, although the bases of this association remain elusive. NK cells have also been implicated in the etiopathogenesis of BD. A family of NK receptors, Killer-cell Immunoglobulin-like Receptor (KIR), with a very complex organization, is very important in the education and control of the NK cells by the union to their ligands, most of them, HLA class I molecules. This study aimed to investigate the contribution of certain KIR functional polymorphisms to the susceptibility to BD. A total of 466 BD patients and 444 healthy individuals were genotyped in HLA class I (A, B, and C). The set of KIR genes and the functional variants of KIR3DL1/DS1 and KIR2DS4 were also determined. Frequency of KIR3DL1*004 was lower in patients than in controls (0.15 vs. 0.20, P = 0.005, Pc = 0.015; OR = 0.70; 95% CI 0.54-0.90) in both B51 positive and negative individuals. KIR3DL1*004, which encodes a misfolded protein, is included in a common telomeric haplotype with only one functional KIR gene, KIR3DL2. Both, KIR3DL1 and KIR3DL2 sense pathogen-associated molecular patterns but they have different capacities to eliminate them. The education of the NK cells depending on the HLA, the balance of KIR3DL1/KIR3DL2 licensed NK cells and the different capacities of these receptors to eliminate pathogens could be involved in the etiopathogenesis of BD.

Diálogos clínicos en el laboratorio: Retos en inmunología del trasplante / No disponible

No disponible

Behçet's disease and genetic interactions between HLA-B*51 and variants in genes of autoinflammatory syndromes.

Behçet's disease (BD) is an immune-mediated systemic disorder with a well-established genetic base. In a previous study, using a next generation sequencing approach, we found many rare variants and some functional polymorphisms in genes related to autoinflammatory syndromes (AID): CECR1, MEFV, MVK, NLRP3, NOD2, PSTPIP1 and TNFRSF1A in our BD cohort. Our strategy did not allow us to establish either number of patients with variants, proportion of individuals accumulating them or relationship with other genetic factors. With the goal to answer these questions, the individual samples were sequenced. Additionally, three functional polymorphisms: NLRP3 p.Gln703Lys, NOD2 p.Arg702Trp and p.Val955Ile were genotyped using TaqMan assays. A total of 98 patients (27.6%) carried at least one rare variant and 13 of them (3.7%) accumulated two or three. Functional regression model analysis suggests epistatic interaction between B51 and MEFV (P = 0.003). A suggestive protective association of the minor allele of NOD2 p.Arg702Trp (P = 0.01) was found in both, B51 positive and negative individuals. Therefore, a high percentage of patients with BD have rare variants in AID genes. Our results suggest that the association of MEFV with BD could be modulated by the HLA molecules; whereas the protective effect of NOD2 p.Arg702Trp would be independent of HLA.

High-resolution characterization of allelic and haplotypic HLA frequency distribution in a Spanish population using high-throughput next-generation sequencing.

Next-generation sequencing (NGS) at the HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1 and -DRB3/4/5 loci was performed on 282 healthy unrelated individuals from different major regions of Spain. High-resolution HLA genotypes defined by full sequencing of class I loci and extended coverage of class II loci were obtained to determine allele frequencies and also to estimate extended haplotype frequencies. HLA alleles were typed at the highest resolution level (4-field level, 4FL); with exception of a minor deviation in HLA-DPA1, no statistically significant deviations from expected Hardy Weinberg Equilibrium (HWE) proportions were observed for all other HLA loci. This study provides new 4FL-allele and -haplotype frequencies estimated for the first time in the Spanish population. Furthermore, our results describe extended haplotypes (including the less frequently typed HLA-DPA1 and HLA-DQA1 loci) and show distinctive haplotype associations found at 4FL-allele definition in this Spanish population study. The distinctive allelic and haplotypic diversity found at the 4FL reveals the high level of heterozygosity and specific haplotypic associations displayed that were not apparent at 2-field level (2FL). Overall, these results may contribute as a useful reference source for future population studies, for HLA-disease association studies as a healthy control group dataset and for improving donor recruitment strategies of bone marrow registries. HLA genotyping data of this Spanish population cohort was also included in the 17th International Histocompatibility and Immunogenetics Workshop (IHIW) as part of the study of HLA diversity in unrelated worldwide populations using NGS.